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We present a systematic method to factorise singularities using non-linear mappings [[!]]. As a 
first application of the method the fully differential NNLO calculation of the partial decay width 
of a Higgs boson into Bottom-quarks is presented [0]. 
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1. Introduction 

With more and more data being collected at the Large Hadron Collider (LHC), it is becoming 
increasingly more important to have control over sizeable higher order QCD corrections to the 
production rates of standard model particles. If theory and experiment are to be compared in any 
meaningful way then these corrections must be computed for arbitrary infrared safe observables 
and cuts on the final state momenta. 

One major difficulty in tackling fully differential calculations beyond the next-to-leading order 
(NLO) is to handle the intricate infra-red divergences in an efficient way. In this talk we show how 
non-linear mappings can be used to factorise such divergences [|]]. 

As a proof of concept we report on the fully differential computation of the H — > bb decay width 
through next-to-next- to-leading-order(NNLO) in the strong coupling expansion ph. 

2. Observables at NNLO in QCD 

For a given observable (O) which may typically be a cross-section or a decay rate, we may 
expand in the strong coupling constant 

(O) = (0 L0 ) + (0 NL0 ) + (^\o NNLO ) + 0(c$) (2.1) 

We will here focus on the NNLO correction (0 NNL0 ), this receives three separate contributions: 

(0 NNW ) = (0 RR ) + (0 RV ) + (O vv ) (2.2) 

Physically the double-real (RR) piece corresponds to the emission of two extra jets in the final 
state, the real-virtual (RV) corresponds to the one-loop corrections as well as one extra jet in the 
final state and the double-virtual (W) corresponds to two-loop corrections. To illustrate these we 
depict interferences of Feynman diagrams corresponding to NNLO corrections to the Higgs to Bot- 
tom quark- anti-quark decay rate in Figure [l| Since the final state phase-space for the RR, RV and 
VV is different, their phase space integrations have to be done separately. The problem which then 
arises is that all three contributions have divergences which only cancel after summing the three 
contributions. These divergences are conveniently regulated in dimensional regularisation, i.e. by 
setting the dimension to d = 4 — 2e, and the problem is thus reduced to finding the Laurent ex- 
pansions of the different NNLO pieces such that after summing them only the finite pieces remain. 
However due to the complicated structure of overlapping divergences at NNLO, this a very non- 
trivial task. Powerful analytical methods exist to deal with loop integrals M, |5[ [7|, 0, |9[ [R]] and 
can also be applied to fully inclusive phase-space integrals. On the contrary dealing with arbitrary 
final state phase-space cuts requires a purely numerical treatment. There are basically two ways in 
which this can be accomplished: 

i) The subtraction method: 

This method is usually based on QCD factorisation to find counter-terms which reproduce the 
limits of amplitudes in the different singular kinematic configurations. Since these counter- 
terms live in kinematic configurations where the concept of infrared-safety does not allow 
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VV: 





RR: 




Figure 1: Sample interferences of Feynman diagrams contributing to the H — > bb decay width 



for experimental cuts to be imposed, they can be integrated analytically and added back to 
the purely virtual piece, such that poles are cancelled analytically. This method has been 



applied successfully in for example [14, H5L 16, 17, h 



ii) The direct integration method: 

This method is more brute force and aims at numerically integrating the real pieces directly. 
One way in which this method has been realized successfully is to use Sector Decomposition 
[11, 12|, 13, 19] to first factorise the singularities and then expand the factorised singularities 
in plus distributions. This may be easily achieved using the following identity: 



-1+ae 



8(x) f, (ae) n 



ae 



rc=0 



log"(jc) 



(2.3) 



where 8 (x) is the Dirac-delta function and the plus-distribution is defined such that 



dxf{x) 



dxg(x) 



/W-/(o) 



Another way to factorise such overlapping singularities is to use Non-linear Mappings [jl|], 
which we will discuss in more detail in the next section. 



3. Factorizing overlapping singularities 

3.1 Sector Decomposition 

We illustrate the method on a simple example integral 



1 x -l+e 



dxdy . (3.1) 

x + y 

To factorise the overlapping singularity at x = = y, we sector decompose by inserting 

l = e(x-y) + 6(y-x) (3.2) 
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where 6 (x) is the Heaviside step function. In each of the two sectors which are created through this 
decomposition we rescale the integration variable back on the unit hypercube, using y t- > xy in the 
first and x t- > xy in the second sector. This factorizes the integrand since both variables contributing 
to the overlapping singularity now scale with either x or y. We obtain 

dxdy- h / dxdy-^- . (3.3) 



lo 1+y Jo l+x 

While this procedure is guaranteed to work for any kind of non-factorised singularity, it will un- 
avoidably lead to a proliferation of integrals. 

3.2 The non-linear mapping 

A single variable mapping which can be used to disentangle overlapping singularities is 

x^P(x,y), P(x,y):= *\ , y>0 (3.4) 

(i x) -\-y 

which keeps the integration boundaries between and 1. This mapping allows, for example, to 
transform the integral I\ into 

r i 

h= dxdy(x y y 1+e (l-x + y y e , (3.5) 



o 



thereby factorizing the singularity at x = = y. It should be noted however that applying y i — >• )3 (y, x) 
yields 

l x~ l+e 

dxdy—— (3.6) 

x+{l-y) 

which does not factorize the singularity but pushes it to the point x = 0,y = 1. This is because the 
denominator of the Jacobian 

dp(x,y) y(l+x) 

dx [(l- x ) +y ]2 V- ' 

is not completely cancelled. Indeed there is a general message to be learned here. First one should 
notice that x already constituted a logarithmic singularity stemming from the factor x~ l+e , while y 
did not since the factor ^ is by itself integrable. 

In [jl|] we therefore termed x to be an active and y a passive singularity. It was then further demon- 
strated that with the successive use of this transformation one can factorize a wide class of dif- 
ferent overlapping structures in multidimensional integrals, including those which typically occur 
in NNLO QCD corrections. In particular it was demonstrated that all double real corrections in 
the hadronic production of arbitrary massive final states and the most singular two-loop integrals 
could be factorised with this mapping. In ^ the mapping was applied to deal with double real and 
real virtual corrections to the process H — > bb. We will discuss some details of this calculation in 
section |j. 

3.3 A two-variable non-linear mapping 

We also mention here a two-variable non-linear transformation 

(x,y) ^ y(x,y), y(x,y) = (xy,^£- ) (3.8) 
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which was suggested to us by Johannes Bliimlein during the RADCOR 201 1 conference. 
Interestingly this mapping can be derived from the following one-variable mapping 

xA 

a(x,A,B), a(x,A,B):=- , (3.9) 

v ' 1 v ' 1 (\-x)B + xA 

which we also employed heavily in [|I|] and [0]. The connection can be seen by noting that 

Y(x,y) = (a(x,a(y,l-x,l),l),a(y,l-x,l)) = (xy *^~^ \ (3.10) 

in other words (x,y) y(x,y) corresponds to the composition of two one- variable mappings, x h-> 
u(x,y, 1) and then y i-> a(y, 1 — x, 1). 

4. NNLO corrections to H -> bb 

All matrix elements which were required for this calculation were generated with QGRAF 
[24]. Further symbolic manipulations, such as color and Dirac algebra, were done with FORM 



[22] and MAPLE [23]. For the computation of loop amplitudes we used the Laporta Algorithm 



[21] implemented in AIR [25] to reduce the amplitude to known master integrals. Numerical 



integrations were done using the VEGAS implementation of the CUBA library [26] 



4.1 Real virtual contribution 

In the real virtual contribution we encounter integrals of type 

2 F 1 (l;-e,l -e,-£) 



d<t>->, 



tu 



where 2F1 (1; — £, 1 — B,x) is Gauss's hypergeometric function, which derives from the all orders in 
£ expansion of the one external mass otherwise massless box. Using a phase-space parametrization 
where the Mandelstam variables, u and t, are factorised it is suggestive to simply expand them with 



the plus-distribution eq.(2.3). However the 2^i(l;— e, 1 — e,x) is singular at x = —00, which is one 



of the subtraction points. This can be seen more clearly when examining the integral representation 



2 F 1 (l,-e,l-e,-l/z)=ze f d y y —^. (4.1) 

Jo z+y 

The integral develops an overlapping singularity when z, — > 0, similar to our earlier example eq.(fO|), 
and can be disentangled with the same mapping. The mapping simply re-derives the well known 
"Pfaff 's transformation" 

2 F l (a,b,c;z) = (l-z)- b 2F l (c-a,b,c;^-). (4.2) 

z—i 

Nevertheless the non-linear mapping should lead to more non-trivial results in cases where such 
identities do not exist. 
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4.2 Double real contribution 

One of the most difficult integrals which one faces in the double real is 

,^ J(Pl,P2,P3,P4) f , A S2a(J(P1,P2,P3,Pa)+J(P3,P4,P1,P2}) a s 
fl<P 4 = / tf<P 4 ■ . (4.3) 

512*34^14^23 J ^34*12*23 (*13*23 + *14*24j 

where J{pi,P2,P3,P4) is a jet function encoding cuts on the final state momenta, and the Lorentz 
invariants are defined such that = Ipi.pj. We use a parametrisation where (we use the shorthand 
Xi = 1 - Xd 



S3A 


= A1A2 




= A1A2A4 


S24 


= A1A2A4 


Sl2 


= A1A2A3 



and 



*13 = h 

514 = Al 



A4A3 + A2A3A4 + 2cos(A57r) y X2X3X3X4X4 



A3A4 + A2A3A4 — 2cos(A57r) y A2A3A3A4A4 



Having partial fractioned and recombined the two terms in the integrand in eq.(|4.3D to avoid the 
line singularity, we end up with a singularity structure similar to 



h 



dxdydz- ' 



(4.4) 



/o "xy(xy + z)' 

The latter can be factorised by first applying x h-> j3(x,z) and then (y,z) i-> (j3(y, 1 — x), j8(z, 1 — x)) 

ii. 

4.3 Numerical results for the H —tbh decay rate 

Our numerical result for the inclusive decay rate is 



r NNLO _ r LO 
1 H^bb 1 H^>bb 



1 + ( -i J 5.6666(4) + ^-i ) 29.14(2) + 0(a*) 



which compares well with the known analytic result [J|] 



r NNLO_ _ -pLO 
1 H^bb L H^bb 



1 + ( 5.6666666.. + (^) 2 29.146714.. + 0(af) 



We also present the 2,3 and 4 jet rates using the Jade algorithm [27] with y cut = 0.01 



r^(4JetRate) = T LO 
//^fc( 3JetRate ) 



--NLO 



•LO 

H^tbb 



+ (^) 2 94.1(l) + ^(a v 3 ) 



+ (^) 19.258(4) + (^) 2 241(2) + G{a]) 



r^(2JetRate) = r"> 



1 



^(2JetRate) - . H ^ bb 
In Figure we also present a fully differential observable. 



f) 13.591(6) -(f) 2 307(2) + ^(a- 



(4.5) 



6 



Disentangling Singularities with Non-linear Mappings 



Franz Herzog 



E max in 2 J et rate with Ycu^ - 1 

lo c:,r:™i 

nlo ==::■.: 

NNLO I I 



0.5 0.51 0.52 0.53 0.54 0.55 

E max /m H 

Figure 2: The maximum energy of the leading jet in the 2 jet rate computed with JADE using y cut =0.1 in 
the H -+bb width. 

5. Conclusion 

In these proceedings we presented a method to factorize singularities at NNLO using a single 
non-linear mapping We showed that the method can be used as a direct integration method 
and thus can be seen as an alternative to Sector decomposition, without proliferating the number 
of integrals. We demonstrated that this method could be used to do entire NNLO calculations by 
presenting the complete fully differential calculation of the H — > bb decay width [Qj. We find that 
our inclusive result is in good agreement with the known analytical result and further present 2,3- 
and 4 jet rates with the JADE algorithm as well as the distribution of the maximum energy of the 
leading jet in the 2-jet rate. 
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